Search for: All records

A new version of the US National Science Foundation National Center forAtmospheric Research (NSF NCAR) thermosphere-ionosphere-electrodynamicsgeneral circulation model (TIEGCM) has been developed and released. Thispaper describes the changes and improvements of the new version 3.0since its last major release (2.0) in 2016. These include: 1) increasingthe model resolution in both the horizontal and vertical dimensions, aswell as the ionospheric dynamo solver; 2) upward extension of the modelupper boundary to enable more accurate simulations of the topsideionosphere and neutral density in the lower exosphere; 3) improvedparameterization for thermal electron heating rate; 4) resolvingtransport of minor species N(2D); 5) treating helium as a major species;6) parameterization for additional physical processes, such as SAPS andelectrojet turbulent heating; 7) including parallel ion drag in theneutral momentum equation; 8) nudging of prognostic fields near thelower boundary from external data; 9) modification to the NO reactionrate and auroral heating rate; 10) outputs of diagnostic analysis termsof the equations; 11) new functionalities enabling model simulations ofcertain recurrent phenomena, such as solar flares and eclipse. Wepresent examples of the model validation during a moderate storm andcompare simulation results by turning on/off new functionalities todemonstrate the related new model capabilities. Furthermore, the modelis upgraded to comply with the new computer software environment at NSFNCAR for easy installation and run setup and with new visualizationtools. Finally, the model limitations and future development plans arediscussed. 

Winds in the nighttime upper thermosphere are often observed to mimic the ionospheric plasma convection at polar latitudes, and whether the same is true for the daytime winds remains unclear. The dayside sector is subject to large temperature gradient set up by solar irradiance and it also contains the cusp, which is a hotspot of Poynting flux and a region with the strongest soft particle precipitation. We examine daytime winds using a Scanning Doppler Imager (SDI) located at the South Pole, and investigate their distribution under steadily positive and negative IMF B_yconditions. The results show that daytime winds exhibit significant differences from the plasma convection. Under negative IMF B_yconditions, winds flow in the same direction as the plasma zonally, but have a meridional component that is strongest in the auroral zone. As a result, winds are more poleward-directed than the plasma convection within the auroral zone, and more westward-directed in the polar cap. Under positive IMF B_yconditions, winds can flow zonally against the plasma in certain regions. For instance, they flow westward in the polar cap despite the eastward plasma convection there, forming a large angle relative to the plasma convection. The results indicate that ion drag may not be the most dominant force for daytime winds. Although the importance of various forcing terms cannot be resolved with the utilized dataset, we speculate that the pressure gradient force in the presence of cusp heating serves as one important contributor. 

Abstract High latitude upper atmospheric inter‐hemispheric asymmetry (IHA) tends to be enhanced during geomagnetic storms, which may be due to the complex spatiotemporal changes and magnitude modifications in field aligned currents (FACs) and particle precipitation (PP). However, the relative contribution of FACs and PP to IHA in high‐latitude forcing and energy is not well understood. The IHA during the 2015 St. Patrick’s Day storm has been investigated using the global ionosphere thermosphere model (GITM), driven by FACs from the Active Magnetosphere and Planetary Electrodynamics Response Experiment (AMPERE) and PP from the Assimilative Mapping of Ionospheric Electrodynamics (AMIE). A comprehensive study of the (a) relative contributions of FACs and PP to electric potential and Joule heating and (b) sensitivity of electric potential and Joule heating to the changes in magnitude and distribution of FACs and PP is presented. The results indicate that FACs lead to larger potential and Joule heating changes compared with PP. The spatial variations of potential and Joule heating are also affected by variation in FACs. As for asymmetric magnitude and distribution, it is found that electric potential and Joule heating are more sensitive to changes in the distribution of FACs and PP than the magnitude of FACs and PP. A new spatial asymmetry index (SAI) is introduced, which reveals spatial asymmetric details that are often overlooked by previous studies. This sensitivity study reveals the relative contributions in high‐latitude forcing and emphasizes the importance of obtaining accurate FACs and PP in both hemispheres. 

The ranked (or top-k) document retrieval problem is defined as follows: preprocess a collection{T₁,T₂,… ,T_d}ofdstrings (called documents) of total lengthninto a data structure, such that for any given query(P,k), wherePis a string (called pattern) of lengthp ≥ 1andk ∈ [1,d]is an integer, the identifiers of thosekdocuments that are most relevant toPcan be reported, ideally in the sorted order of their relevance. The seminal work by Hon et al. [FOCS 2009 and Journal of the ACM 2014] presented anO(n)-space (in words) data structure withO(p+klogk)query time. The query time was later improved toO(p+k)[SODA 2012] and further toO(p/log_σn+k)[SIAM Journal on Computing 2017] by Navarro and Nekrich, whereσis the alphabet size. We revisit this problem in the external memory model and present three data structures. The first one takesO(n)-space and answer queries inO(p/B+ log_Bn + k/B+log^*(n/B)) I/Os, whereBis the block size. The second one takesO(nlog^*(n/B)) space and answer queries in optimalO(p/B+ log_Bn + k/B)I/Os. In both cases, the answers are reported in the unsorted order of relevance. To handle sorted top-kdocument retrieval, we present anO(nlog(d/B))space data structure with optimal query cost. 

Inter-hemispheric asymmetry (IHA) in Earth’s ionosphere–thermosphere (IT) system can be associated with high-latitude forcing that intensifies during storm time, e.g., ion convection, auroral electron precipitation, and energy deposition, but a comprehensive understanding of the pathways that generate IHA in the IT is lacking. Numerical simulations can help address this issue, but accurate specification of high-latitude forcing is needed. In this study, we utilize the Active Magnetosphere and Planetary Electrodynamics Response Experiment-revised fieldaligned currents (FACs) to specify the high-latitude electric potential in the Global Ionosphere and Thermosphere Model (GITM) during the October 8–9, 2012, storm. Our result illustrates the advantages of the FAC-driven technique in capturing high-latitude ion drift, ion convection equatorial boundary, and the storm-time neutral density response observed by satellite. First, it is found that the cross-polar-cap potential, hemispheric power, and ion convection distribution can be highly asymmetric between two hemispheres with a clear B_ydependence in the convection equatorial boundary. Comparison with simulation based on mirror precipitation suggests that the convection distribution is more sensitive to FAC, while its intensity also depends on the ionospheric conductance-related precipitation. Second, the IHA in the neutral density response closely follows the IHA in the total Joule heating dissipation with a time delay. Stronger Joule heating deposited associated with greater high-latitude electric potential in the southern hemisphere during the focus period generates more neutral density as well, which provides some evidences that the high-latitude forcing could become the dominant factor to IHAs in the thermosphere when near the equinox. Our study improves the understanding of storm-time IHA in high-latitude forcing and the IT system. 

An important question that is being increasingly studied across subdisciplines of Heliophysics is “how do mesoscale phenomena contribute to the global response of the system?” This review paper focuses on this question within two specific but interlinked regions in Near-Earth space: the magnetotail’s transition region to the inner magnetosphere and the ionosphere. There is a concerted effort within the Geospace Environment Modeling (GEM) community to understand the degree to which mesoscale transport in the magnetotail contributes to the global dynamics of magnetic flux transport and dipolarization, particle transport and injections contributing to the storm-time ring current development, and the substorm current wedge. Because the magnetosphere-ionosphere is a tightly coupled system, it is also important to understand how mesoscale transport in the magnetotail impacts auroral precipitation and the global ionospheric system response. Groups within the Coupling, Energetics and Dynamics of Atmospheric Regions Program (CEDAR) community have also been studying how the ionosphere-thermosphere responds to these mesoscale drivers. These specific open questions are part of a larger need to better characterize and quantify mesoscale “messengers” or “conduits” of information—magnetic flux, particle flux, current, and energy—which are key to understanding the global system. After reviewing recent progress and open questions, we suggest datasets that, if developed in the future, will help answer these questions. 

Abstract Meso‐scale plasma convection and particle precipitation could be significant momentum and energy sources for the ionosphere‐thermosphere (I‐T) system. Following our previous work on the I‐T response to a typical midnight flow burst, flow bursts with different characteristics (lifetime, size, and speed) have been examined systematically with Global Ionosphere‐Thermosphere Model (GITM) simulations in this study. Differences between simulations with and without additional flow bursts are used to illustrate the impact of flow bursts on the I‐T system. The neutral density perturbation due to a flow burst increases with the lifetime, size, and flow speed of the flow burst. It was found that the neutral density perturbation is most sensitive to the size of a flow burst, increasing from ∼0.3% to ∼1.3% when the size changes from 80 to 200 km. A westward‐eastward asymmetry has been identified in neutral density, wind, and temperature perturbations, which may be due to the changing of the forcing location in geographic coordinates and the asymmetrical background state of the I‐T system. In addition to midnight flow bursts, simulations with flow bursts centered at noon, dawn, and dusk have also been carried out. A flow burst centered at noon (12.0 Local Time [LT], 73°N) produces the weakest perturbation, and a flow burst centered at dusk (18.0 LT, 71°N) produces the strongest. Single‐cell and two‐cell flow bursts induce very similar neutral density perturbation patterns. 

Abstract In this study, the Global Ionosphere Thermosphere Model is utilized to investigate the inter‐hemispheric asymmetry in the ionosphere‐thermosphere (I‐T) system at mid‐ and high‐latitudes (|geographic latitude| > 45°) associated with inter‐hemispheric differences in (a) the solar irradiance, (b) geomagnetic field, and (c) magnetospheric forcing under moderate geomagnetic conditions. Specifically, we have quantified the relative significance of the above three causes to the inter‐hemispheric asymmetries in the spatially weighted averaged E‐region electron density, F‐region neutral mass density, and horizontal neutral wind along with the hemispheric‐integrated Joule heating. Further, an asymmetry index defined as the percentage differences of these four quantities between the northern and southern hemispheres (|geographic latitude| > 45°) was calculated. It is found that: (a) The difference of the solar extreme ulutraviolet (EUV) irradiance plays a dominant role in causing inter‐hemispheric asymmetries in the four examined I‐T quantities. Typically, the asymmetry index for the E‐region electron density and integrated Joule heating at solstices with F10.7 = 150 sfu can reach 92.97% and 38.25%, respectively. (b) The asymmetric geomagnetic field can result in a strong daily variation of inter‐hemispheric asymmetries in the F‐region neutral wind and hemispheric‐integrated Joule heating over geographic coordinates. Their amplitude of asymmetry indices can be as large as 20.81% and 42.52%, which can be comparable to the solar EUV irradiance effect. (c) The contributions of the asymmetric magnetospheric forcing, including particle precipitation and ion convection pattern, can cause the asymmetry of integrated Joule heating as significant as 28.43% and 34.72%, respectively, which can be even stronger than other causes when the geomagnetic activity is intense.